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Abstract: Spin-2 fields are often candidates in physics beyond the Standard Model 
namely the models with extra-dimensions where spin-2 Kaluza-Klein gravitons couple to 
the fields of the SM. Also, in the context of Higgs searches, spin-2 fields have been studied 
as an alternative to the scalar Higgs boson. In this article, we present the complete three 
loop QCD radiative corrections to the spin-2 quark-antiquark and spin-2 gluon-gluon form 
factors in SU(N) gauge theory with nj light flavors. These form factors contribute to 
both quark-antiquark and gluon-gluon initiated processes involving spin-2 particle in the 
hadronic reactions at the LHC. We have studied the structure of infrared singularities in 
these form factors up to three loop level using Sudakov integro-differential equation and 
found that the anomalous dimensions originating from soft and collinear regions of the 
loop integrals coincide with those of the electroweak vector boson and Higgs form factors 
confirming the universality of the infrared singularities in QCD amplitudes. 
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1 Introduction 

In the context of the recent discovery of the new boson at the LHC, with mass of about 125 
GeV [1, 2], there has been renewed interest in massive spin-2 resonance which could also 
lead to similar final states [3]. The massive spin-2 could be a Kaluza-Klein (KK) graviton 
of the TeV scale gravity models [4, 5] as a result of gravity propagating in the extra 
dimensional bulk or any generic spin-2 resonance in some other new physics scenarios. It 
was noted in [6] that gauge symmetry and Lorentz invariance forbid operators of dimension 
four that could lead to a coupling of a massive spin-2 resonance to a pair of SM particles. 
Further, if the flavor and CP symmetries of the SM are respected by these new physics 
scenarios, the leading dimension five operator is none other than the energy momentum 
tensor of the SM particles. The structure of the operator coupling thus being identical 
to the KK graviton, though the constant coefficients could be different for the KK graviton 
or any generic spin-2 imposter. Nonetheless, methods to distinguish KK graviton from the 
imposter have been proposed [6] and will be of importance for BSM searches at the LHC 
which is now operational at higher energies and luminosity. The increasing accuracy of the 
experimental data at the LHC Run-II, demands an equally precise theoretical predictions. 

To match the current theoretical accuracy of say the Drell-Yan production [7-9], the 
Higgs boson production in gluon fusion [10-16], in bottom quark annihilation [17, 18] 
and associated production with a vector boson [19, 20] at the LHC, it is imperative that 
competing BSM models are also available to the same accuracy in higher orders in QCD. 
Form factors are essential ingredients for many precision calculation in QCD. An important 
building block for phenomenological study is the computation of qq —» spin-2 and gg —)> 
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spin-2 form factors and is at present available to two-loop in QCD [21], while for many 
processes of interest they are now available to the 3-loop order [22-27]. Stringent bounds 
[28-30] on the parameters of ADD and RS models are available due to the presence of 
precise theoretical predictions for various important observables up to NLO level in QCD. 
Often, these observables suffer from large uncertainties resulting from renormalization and 
factorization scales and the only remedy to this is to include higher order QCD effects 
to the born contributions. The NLO QCD predictions based on fixed order as well as 
parton shower improved in the MadGraph5_aMC@NLO [31] framework for di-final state 
[32-40] productions in the gravity mediated models have already played important role in 
constraining the model parameters of ADD and RS. Next-to-next-to leading order (NNLO) 
corrections for the graviton production at the LHC in the threshold limit are already 
available [41] and attempts to improve these predictions through NNLO corrections and 
beyond are already underway [42]. As these corrections are only sensitive to the tensorial 
interaction and not sensitive to the details of the model, these results are applicable to 
production of any generic spin-2 resonance. Hence this article takes the first step towards 
going beyond NNLO for the resonant production of a generic spin-2 particle at the LHC, 
namely the computation of quark and gluon form factors at three loop level in perturbative 
QCD with rtf light flavors. We will report the first results on the threshold effects at N 3 LO 
in the future publication and demonstrate the importance such corrections at the LHC in 
the context of spin-2 resonance searches. 

In addition to the phenomenological importance with respect to precise predictions 
of some observable, form factors in QCD are of considerable theoretical interest in terms 
of the factorization and universal nature of the singular structure. Studying the infrared 
pole structure and factorization properties of these IR singularities in multi-loop QCD 
amplitudes with tensorial coupling to 3-loop order and to confirm the standard expectation 
of QCD amplitudes [43-46] is an essential prerequisite. The spin-2 field being a tensor of 
rank-2 is coupled to the energy-momentum tensor T^, which is a symmetric and conserved 
quantity. The operator of QCD is finite [47], which would imply no UV renormalization 
is required. Further consists of gauge invariant terms and in addition gauge dependent 
and ghost terms, we explicitly observe that to the 3-loop order these form factors are 
independent of the gauge dependent and ghost terms [47, 48], which is an important check 
of the calculation. From a computational point of view 3-loop amplitudes with higher 
tensorial coupling is being attempted for the first time. At the intermediate stages of 
the computation this leads to higher rank tensorial integrals resulting from more than 
3000 three loop Feynman amplitudes contributing to the gluon form factor alone. This 
computation again establishes the power of several state-of-the-art techniques namely IBP 
and LI identities and differential equation method to solve the master integrals. 

In the next section 2, we describe the effective Lagrangian. In section 3, after defining 
the quark and gluon form factors, we present the computational details at three loop level 
followed by the results. The details of ultraviolet renormalization and universal structure 
of infrared poles are given in section 4 and section 5 respectively. Finally we conclude with 
our findings in section 6. 
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2 The Effective Lagrangian 


The effective Lagrangian that describes the interaction of the spin-2 field with the SM 
fields can be written down in a gauge invariant way through the energy momentum tensor 
of the SM fields. We denote the spin-2 field by h^ v and the SM energy momentum tensor 
by . Since we are interested only in the QCD corrections to processes involving spin-2 
fields, we restrict ourselves to the QCD part of and the corresponding action reads 
[4, 5] as 


5 = 


Ssm + s h - ^ J d A x T®° d (x) h txu {x ), 


(2.1) 


where Ssm is the SM action, Sh is the kinetic energy part of the action corresponding to 
spin-2 fields, k, is a dimensionful coupling and T^ D is the energy momentum tensor of 
QCD given by 

Tfr CD = -g^CQCD - f; p F^ - -^ p dP{A a p d"Al) + + A°d v {6rA%)) 

- ig s T a Al)ip - 'ipCdu + ig s T a Al)^ p il> + Vnv( 7^ - igsT a A a p )i) 
-?(K + ig s T a A“) 7^1 + d p oJ a (d„uj a - g s f abc A c u u b ) 


l 

+ 4 


+d u uf(d p u; a - g s f abc A c p J). 


( 2 . 2 ) 


g s is the strong coupling constant and £ is the gauge fixing parameter. The T a are genera¬ 
tors and f abc are the structure constants of £77(3). Note that spin-2 fields couple to ghost 
fields (c j a ) [49] as well in order to cancel unphysical degrees of freedom of gluon fields (A^). 


3 The Form Factors 

The form factor parametrizes the interaction of the spin-2 field with those of SM order by 
order in perturbation theory. We compute both quark and gluon form factors by sand¬ 
wiching the energy-momentum tensor between on-shell quark and gluon states respectively 
normalized by their respective born amplitudes: 


Fl(Q\e) = 


.m| 0) *..m1 0) 


l I 

oo 

71—0 


s 


5 e n jf’ (n) (e), 


I = q,g 


(3.1) 


where M\ are the unrenormalized amplitudes computed in powers of the bare strong cou¬ 
pling constant a s = g^/lftn 2 using dimensional regularization in d = 4 + e dimensions, that 





where Q 2 = — 2p\ • p 2 and pi, P 2 are the momenta of external quark or gluon on-shell states. 
The dimensionful scale p is introduced to keep the strong coupling constant dimensionless 
in d space-time dimensions. The other constant at n th loop is S™ = exp[^( 7 £ — ln 47 r)] 
where Euler constant 7 e = 0.5772.... 

In [21], both one and two loop form factors were presented in dimensional regularization 
and later on, they were used in [41] to compute the threshold corrections to Drell-Yan 
production at the LHC in ADD and RS models to second order in strong coupling constant. 
In the following, we present the third order correction to the form factors in QCD. 


3.1 Computational Procedure 

In this section, we describe in detail the method that we follow to compute both quark 
and gluon form factors of the energy momentum tensor to third order in strong coupling 
constant using dimensional regularization. The relevant amplitudes are generated using 
QGRAF [50]. At third order alone, there are 3374 and 1072 number of Feynman diagrams 
for gluon and quark form factors respectively. The QGRAF generated amplitudes are then 
converted into a suitable format using routines developed using the symbolic manipulation 
program FORM [51]. Both group as well as Lorentz indices are carefully handled to 
express the form factors in a suitable color basis involving Casimir operators of SU(N ) 
with the coefficients containing three loop scalar integrals. For the gluon form factor we 
have summed only the physical polarizations of the external gluons using 


^e^(pi,s)e u *(pi,s) 


-<r + 


gM+ qM 

Pi-Qi 


(3.3) 


where, pi is the i^'-gluon momentum and qi is the corresponding light-like momentum. We 
choose qi = p‘> and </■_> = Pi for simplicity. For the external spin-2 fields, we have used the 
d dimensional polarization sum given in [32] with q being the spin-2 momentum 


B^ p<T (q) 



ql'qP\ 
<M / 



<fq°\ 
<M / 


+ 



qP (f \ 
<M / 



qW_ \ 

<M ) 


2 

d — 1 



<fq*\ 
q-q ) 



qpqp\ 

q-q ) 


(3.4) 


We have used Feynman gauge throughout. 

At three loop level, we find that the diagrams contributing to form factors can have at 
most 9 independent propagators involving two external momenta p \, p 2 and three internal 
loop momenta k\,k 2 ,k^, while the maximum number of scalar products that can appear 
in the numerator of each diagram can be 12. Hence we need to increase the number of 
propagators to 12 which allow us to classify all the three loop diagrams into three different 
auxiliary topologies. We take the help of Reduze2 [52] for this purpose. The topologies 
[25] that are used in our computation are given below 


Ai : 1*1,1*2, T , 3,T , i2,T ) 13,T , 23,^l;l)^l:12,T>2;l, 1*2:12, T , 3 ; i,P3:12 

A2 : Pi, P 2 , Z>3, Z>12, Z>13, ^23, X>13;2, ©1:12, ©2;1, ®12:2, Z>3;1, ©3:12 
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A 3 : T>i,V2,V3,Vi2i'Di2,'Di23i'Di- ) i,Vi : i2,'D2,l,'D2:12i'D3 ) li'D3:i2 


(3.5) 


where, 


Vi = k 2 , V i:j = (ki - kj) 2 , V i;j i = (.ki - kj - k t ) 2 , 

T>i-,j = (h ~ Pj) 2 , Vi-jl = (h - pj - pi) 2 , Vij-i = (ki - kj - pi ) 2 . (3.6) 


The resulting integrals classified in terms of three topologies, are then reduced to a set 
of master integrals by using a systematic approach that uses Integration by parts (IBP) 
[53] and Lorentz invariant (LI) [54] identities. The IBP identities follow from the fact 
that within dimensional regularization, the integrals are finite and well-behaved and hence 
any integrand at the boundary must be zero. Following this, the generalization of Gauss 
theorem implies the integral of the total derivative with respect to any loop momenta to 
be zero, that is 


/ 


d d k\ 

{2n) d 


d d ks d 
(! 2n) d dki 



(3.7) 


where ni is an element of n = (ni, • • -, 7112 ) with rq E Z and Dis are propagators which 
depend on the loop and external momenta. The four vector Vj can be both loop and 
external momenta. Performing the differentiation on the left hand side and expressing the 
scalar products of ki and pj linearly in terms of V{ s, one obtains the IBP identities given 
by 

]>>J(^i + ni,...A,i2 + ni2) = 0 (3.8) 

i 


where 

J(m) = Amu--;m l2 ) = J ^ (3.9) 

with bij G { — 1, 0,1} and ai are polynomial in rij. The LI identities follow from the fact that 
the loop integrals are invariant under Lorentz transformations of the external momenta, 
that is 

P?Pj J (”) = °- ( 3 - 10 ) 

V k vPk / 

For the case of three loop form factor, there are 15 IBP identities and 1 LI identity for 
each integrand, and hence there are large number of equations for the whole system. These 
equations can be solved to relate the large number of scalar integrals and express them 
in terms of a set of fewer integrals which are the so called master integrals. To solve this 
large system of equations, there are dedicated computer algebra tools like AIR [55], FIRE 
[56], REDUZE [52, 57], LiteRed [58, 59] etc. We use the Mathematica based package 
LiteRedV1.82 along with MintVI. 1 [60]. 

We find that the form factors at three loop level can be expressed in terms of 22 master 
integrals. Following the same notation as of [25], the master integrals can be distinguished 
into three topological types: genuine three loop integrals with vertex functions (A^), three 
loop propagator integrals (B t ^) and integrals which are product of one loop and two loop 
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integrals (C^). Defining a generic three loop master integral as 


A- 



i = 1,2,3 (3.11) 


where Dj is the j th element of the set Ai, we can identify the resulting master integrals in 
our computation with those given in [25] and they are listed 1 in Fig. 1 and Fig. 2. 





£?6,1 = ^1,111000010101 


^>6,2 = ^1,011110000101 £>8,1 = ^3,011111010101 



c 6tl = ^1,011100100101 Cg,l = ^2,111100011101 

Figure 1: Two point and factorizable three point three loop integrals 


The master integrals were computed in [24, 25] to relevant orders in e and we have 
used them to complete our computation of the form factors up to three loop level. The 
electronic version of the results of both quark and gluon form factors in terms of the master 
integrals Aij^Bij and Cij for arbitrary d is attached with the arXiv version. In the next 
section, we present the three loop results for both the form factors expanded in powers of 
e along with already known one and two loop results. 


3.2 Results 

In this section we present one, two and three loop quark and gluon form factors after 
expanding in powers of e to relevant order. The one and two loop results completely agree 
with [21] and the three loop ones are new. 


fT,( l) 

^9 



203 

TsT 
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37307 77„ 

2592 36 s 


1 These flgure have been taken from [25]. 


- 6 - 










^- 5,1 = ^- 1,001101100001 ^- 5,2 = ^- 1,001011011000 ^- 6,1 = ^- 1,010101100110 












Figure 2: Genuine three loop integrals with vertex function 
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where Ca = iV, Cp = ( N 2 — 1)/2N and rif is the number of active quark flavors. In the 
next section, we describe how these form factors can be renormalized up to three loop level 
through coupling constant renormalization. We then study the universal structure of the 
infrared poles in e through Sudakov’s KG equation up to three loop level. It provides a 
crucial check for our new results on the form factors. 


(3.16) 


4 Ultraviolet Renormalization 


In MS scheme, the renormalized coupling constant a s = a s (ii%) at the renormalization 
scale hr is related to unrenormalized coupling constant a s by 


= “T t 1 + a * r i + 2 + C’(af)] 

M R 


(4.1) 


where, 


2 A) / 4/3q Pi 
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where Pi are the coefficients of QCD beta function : 


A) = ( yCU - , /3i = (jC 2 a - ™C A n f - 2 C F n f 


(4.2) 


Using the Eq. 4.1, we now can express M\ (Eq. 3.1) in powers of renormalized a s with UV 

( i ) 

finite matrix elements A4j 
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Using above equations, we can obtain the renormalized form factors J -j T in terms of a s . 
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5 Infrared Singularities and Universal Pole Structure 


The results on multiparton amplitudes beyond leading order in perturbative QCD have not 
only played an important role in understanding the infrared structure of the theory but also 
allowed us to successfully carry out various resummation programs for physical observables 
in the kinematic regions where the fixed order perturbation theory breaks down. The most 
important one along this line was the very successful proposal by Catani [43] (see also 
[44]) on one and two loop QCD amplitudes using the universal subtraction operators. The 
generalization of this proposal was achieved by Becher and Neubert [45] and also by Gardi 
and Magnea [46] beyond two loops. In [61], for the first time, the structure of single pole 
term in both quark and gluon form factors at two loop level was unraveled. It was shown 
explicitly that the single pole can be written as a linear combination of UV, collinear and 
soft anomalous dimensions. The fact that this feature continues to hold even at three loop 
level for the same form factors was observed in [22]. The structure of the single pole term 
for the multiparton amplitudes was studied in detail in [62, 63]. 

The form factors Jj (a s , Q 2 , /i 2 , e) satisfy the KG-differential equation which follows 
from the factorization property, gauge and renormalization group invariances [64-67] 


Q 2 (a s , Q 2 , 1 ?, e) = I 


r/T.I/;, Hr , r T,I/» Q 2 I j2 R c \ 
M 2 Hr H 2 


(5.1) 


where, all the poles in dimensional regularization parameter e are contained in A- which 
is also taken to be Q 2 independent and the finite terms as e —> 0 are encapsulated in G-. 
Renormalization group invariance of the form factor implies 


Hr 




K 


T,I 


(a s , ^§,e) = 
H 


-Hr 


dH 2 R 


r T ’ l (n 9l VR A = 
^ lAsi O’ 2 1 G 

Hr H 2 


■J2 a s(HR)A 2 


T,I 


(5.2) 


2=1 


where, A T are the cusp anomalous dimensions. Since K 1 in Eq. 5.2 contains only poles in 
e with no Q 2 dependence, it can be easily solved in powers of bare strong coupling constant 
a s . Expressing 


K 


T.I 




H 


i = 1 


Hr 


H 


(5.3) 


we find that the constants K T 4;W( e ) consist of simple poles in e with the coefficients 
containing Aj jI ’ and s. These can readily be found in [68, 69]. 

The renormalization group equation of G T,I (a s , ,e) can be solved by imposing 

Hr H 

the boundary condition at /jL 2 r — Q 2 . Hence the solution can be expressed in terms of the 
boundary function G T,I (a s (Q 2 ), l,e) and the term that contains full dependence : 


G T ’ 1 (a s , , e) = G T,1 (a s (Q 2 ] 

Hr H 


, 1 , e) + 


[ Q 2 ^r^ T,I ( a s(^l)) (5-4) 

/ w X 

d —o- 
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(5.5) 


The fi 2 R independent part of the solution can be expanded in powers of a s as 

oo 

G T ’\a s (Q 2 ), 1, e) = f] ai(Q 2 )Gj’\e ). 

2=1 

Substituting Eq. 5.3 and Eq. 5.4 in Eq. 5.1, and integrating over Q 2 , we obtain the form 
factor in powers of strong coupling constant: 


00 /n 2 \ i 2 

In I-[(a s ,Q 2 ,fi 2 ,e) = J2^[%) %£}$> (e) 

i = 1 W / 


(5.6) 


with 




-r r 

* 1 ,( 2 ) 


1 | Til II 1 t I n syT I / \ I 1 | 1 


eg’if) = irA /%A, ,J [tjj- 5A" - AGf’‘( £ ) f ^ 


€ 2 


4 ?^) = 1 +1 '/Ml 1 + ? 


+ - i^l 1 - jft at'(e) - V + ^ 2G‘-‘( 


VT,I 


T,I / 


1 1 


vr.i / 


9 


e 3 


(5.7) 


It is now straightforward to extract the cusp anomalous dimensions by comparing Eq. 5.7 
with the form factors presented in the previous section. We find 


AV=C,{l). 


A2’ 1 = CiCa \ - 8C2| + Cirif • 

2 , 490 1072C2 88C3 176C| 1 


^ = CiC\\^f 


, 88 C 3 , 176C 2 2 1 , ^ „ f HO , 00 ^ 1 
+ ^- + — j + OCVn/ 3- + 32Cs| 


„ „ [ 836 I 6 OC 2 112Cs 1 ~ 2 f 16 

+ CiC A n, - -=} + C,nj 

where C\ — Cp for I — q and C\ — Ca for I — g. We find that they not only satisfy the 
property of maximally non-abelian but also coincide with those that appear in the quark 
and gluon form factors which are available up to three-loop level in the literature [70, 71]. 


A J ’ 9 = 7 T iA J’ q and A J ’ 1 = A i I = Q , 9 ■ 
Of 

T I 

Following [22] and [61], we can parametrize G i ’ (e) as follows: 


(5.8) 


OO 


et T.m 


(5.9) 


k =1 
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T I 

where the constants C i ’ are given by [69] 


cJ l = o, 

cf' = -20 o9 ^m 

cf'- 


= -2A ! ,^' 1 »-2A,( 9 J« 1) +2A 9 y«) 


(5.10) 


Using the above decomposition, we can extract dJ’ 1 and fj' 1 from the form factors com¬ 
puted up to three loop level. They are found to be 


Bj* = C A 


~ n fU 


Bl 9 = C 2 a{j + 12Cs} - n f C A ||| - n f C F { 2 } , 
bJ’ 9 = C A C F n f 


B^ q = 


2411 ^ 2 [291 2 f 233 8 A 4 a2 80 . 1 

L ^ n/ \l^ + 3 C2 + 3 C2 + T C3 } 

+ C a{y~ + |C2 + f Cl + ~ 80Cs} + C F nj {y } + C 2 F n f { l} , 


Cf|i| 1 


Bl' 1 = Cf\1- 12 C2 + 24(3 \ + 


Bj’ q = C\C F { - 2 C 2 2 + ^C2 - + 40Cs 


C aC F \jL + ? C2 12Cs } + j -1 - y <C2 j , 


9 


1657 

~36~ 


^C2 + fCs + 120(5 + 1 ~f\+ C A C F n f \ ^2 2 


^ 1 

988 . 

+ 


1336 A 

200 A , 

27 C2 + 

X“ Cs + ^ 


| ~ 32 C 2 C 3 + I8C2 + 68C3 - 240 C 5 + "2" | 

„ 2 ( 232 2 20 136 1 ^ 2 f 80^ 16 A 171 

+ 2 + Y C2 - ^-C3 - 231 + C F n f ^-C 2 - yCs - y j • 


(5.11) 

that appear in quark and gluon form 


x 1 

We find that the above B { ’ are identical to the ones 
factors of [22]: 

Bj J = B \, / = 9 , 5 , i = 1,2,3 (5.12) 


and 


/i T,I = o, 


/ 2 tj 


f T,l (-".52 , 170 . 12650 1316 136781 ] 

h =C t C A |—C2 + x C2C3-^ i -< 2 -^-<3 + 192C 5 +^-j 
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+ CiC A n f 


11842 

729 


96 . o 2828 A 728 A 

“T C2 + ^T C2 + ^7 C3 “' 

,32 o _ 304 17111 

+ ClC F n f \ —( 2 + 4C2 + -g-C3- 

„ ZJ 40 112 20801 

+ i n f \ 27^ + 27 ^ 729 } ' 


(5.13) 


Similar to cusp anomalous dimensions, we find that fj' 1 satisfy the property of maximally 
non-abelian and in addition they coincide with those that appear in the quark and gluon 
form factors which are available up to three-loop level in the literature [61], 


/7 s = ^/7’ /7' = fi 


I = q ,9 • = 1,2,3. (5.14) 


The UV anomalous dimensions are found to be identically zero due the conservation of 
QCD energy momentum tensor, i.e., 


t,i n 

V = o. 


(5.15) 


The universal behavior of infrared poles in terms of the cusp (A 1 ), collinear (B 1 ) and soft 
(/*) anomalous dimensions provides a crucial check on our computation. The remaining 
terms namely g j n gq 5.9 can be extracted from the form factors and they are listed 

below: 




T r( 2 ) ^ / 2879 7 1L \ / 497 lW 

9l ~ Ca [ Vl 6 ~ ~~ 3^ ~~ 12 ^ ) + n/ ( ~~ 216 + 6 ^ ) ’ 


9 1 

92 


T,g ; (3) ^ ( 37307 77 203 47 A ^6593 7 A 35 , 

Ca ( 2592 + 36^ + 144^ 2 + 80) + Hf ^2592 18^ 144 ^ ) ’ 


T ’ s;(1) = Cl - 


19333 88 . 799 „ \ ^ (34991 32 82 P 

~ur + Y Cs + T^ C2 +CAn n^2T + Y Cs - ¥ C2 


61 


+ C F n f \ — - I 6 C 3 ) +n f | 


2219 2 \ 

~324~ + 9^ ) ’ 


92 


3888 


448 183 221 P\ _ ( 2245 118 . . 24 A 

"¥ (3 + T 2 '¥ (2 j+CpnA -^- + ^-C3 + C 2 + yC 2 


2 /1999 14 35 \ 

+ n f I T77T - ^Cs - ^C 2 ) , 


162 27 


36 


(5.16) 


15 - 

















!/,“'' ' 1 ' = C F I - 10 + < 2 ) , 


£" m =C r (l2-| 6 -| & ). 


T,q;( 3 ) _ n ( 1Q ,7 5 47 

9l “ 13+ 4 C3 + 4 C2 + 80 C2 ’ 


#2 


T,q;(i)_<o2 ( 107 10yl/ , , nn , 88 a2 \ ~ ~ / 91693 . 452 . 1103 


“T2"“ 124C3 + 90C2 “T C2 ) +CaCf [ “^r + ^“ C3 “ 18 


C 2 


92 


88 2 \ / 7397 8 85 \ 

+ y C2 J + C F n f - -Cs + y C2 J , 

T, q;( 2) = cl ( 12 ^ + 12C5 + 328Cs - ^C2 - 28C2C3 + ^Ci 


V 


48 


12 


15 


(2192269 riA 20399 . 15751 89 . . 2027 a2 

+ A F y 3888 ^ 54~^ 3 + 108 + T C2Cs ~~ 40 ^ 

„ ( 168557 59. 1079. 7 >2 \ 

+ CFn A~^M T-27 C3 -^T C2+ 12 C2 ' 


(5.17) 


6 Conclusions 

We have presented both quark-antiquark and gluon-gluon form factors of the spin-2 fields 
that couple to fields of SU(N) gauge theory with rtf light flavors. We have used state- 
of-the-art methods to perform this computation efficiently as the number of Feynman 
diagrams involved is quite large compared to other known form factors. We have used 
IBP and LI identities to express the form factors in terms of 22 master integrals. We have 
presented the form factors in terms of these master integrals for arbitrary d as well as in 
powers of c — d — 4 to appropriate order, thanks to the availability of the master integrals 
to relevant orders in e for further study. These form factors are important components 
to the scattering cross sections involving spin-2 fields beyond leading order in QCD. We 
have shown that these form factors do satisfy Sudakov integro-differential equation and 
hence exhibit identical infrared structure of other form factors such as those appearing 
in electroweak vector boson and Higgs productions up to three loop level. We have also 
shown these factors do not require overall renormalization due to the conservation property 
of the energy momentum tensor. Our results will be useful in improving the perturbative 
predictions of spin-2 resonance production beyond NNLO level at the LHC where searches 
for such particles are already underway with the upgraded energy and luminosity. 
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